Ground-State Structure in z/ = 2 Bilayer Quantum Hall Systems 
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We investigate the ground-state structure of the bilayer quantum Hall system at the filling factor 
V = 2. Making an exact analysis of the ground state in the SU(4)-invariant limit, we include all other 
interactions as small perturbation. We carry out analytic calculations and construct phase diagrams 
for nonzero values of the Zeeman, tunneling and bias interactions. In particular we examine carefully 
how the phase transition occurs by applying the bias voltage and inducing a density imbalance 
between the two layers. We compare our theoretical result with the experimental data due to 
Sawada et al. based on the phase diagram in the (tq-Po plane, where po and cro are the total electron 
density and the density difference between the two layers, respectively. 



I. INTRODUCTION 

A rich physics has emerged from the layer degree of 
freedom in the bilayer quantum Hall (QH) systemii^. 
An electron carries the SU(2) pseudospin index assigned 
to the front and back layers. The spin and the pseu- 
dospin are equal partners. There arises a unique phase 
at the filling factor v — \. It is the spin-ferromagnet 
and pseudospin- ferromagnet phase, where various in- 
tralayer and interlayer coherent phenomena have been 
observedii^. On the other hand, according to the one- 
body picture we expect to have two phases at v = 2 
depending on the relative strength between the Zeeman 
gap Az and the tunneling gap Asas- One is the spin- 
ferromagnet and pseudospin-singlet phase (abridged as 
the spin phase'^) for Az > AgAs; the other is the spin- 
singlet and pseudospin-ferromagnet phase (abridged as 
the ppin phase'^) for Az < Asas- Between these two 
phases, driven by interlayer correlations, a novel canted 
antiferromagnetic phase has been predicted to emerge^. 

The first experimental indication of a phase tran- 
sition at = 2 was revealed by inelastic light 
scattering spectroscopj»&. An unambiguous evidence 
was subsequently obtained through magneto-transport 
measurements^, where the importance of making density 
imbalance between the two layers is emphasized to study 
phase transitions. There is an attempt? to interpret the 
data— based on numerical analyses. However, the phys- 
ical picture is not yet clear enough due to the lack of 
analytic understanding. For instance, it is yet to be ex- 
plored how the density imbalance is induced as a function 
of the bias voltage Fbias- The density imbalance is made 
as soon as Vbias ^ Q aX v = 1, but this is not so simple at 
v = 2. Furthermore, the exact diagonalizationS^ of a few- 
electron system has shown that the boundary between 
the spin and canted phases is practically unmodified but 
the boundary between the canted and ppin phases is con- 
siderably modified from the mean- field result. There is 
so far no theoretical explanation of these behaviors. The 
aim of this paper is to understand the ground-state struc- 
ture of the V = 2 bilayer QH system more in details based 
on perturbation theory. 

We investigate physics taking place in the lowest Lan- 



dau level (LLL). Since each Landau site can accommo- 
date four electrons with the spin and pseudospin degrees 
of freedom, the underlying group structure is enlarged to 
SU(4). We start with the SU(4)-invariant hmit of the 
bilayer system, where the ground state is determined by 
solving the eigenvalue equation. We then include SU(4)- 
noninvariant terms perturbatively. The ground-state en- 
ergy is calculated analytically in the first order of per- 
turbation. Then we examine carefully how the phase 
transition occurs as the density imbalance is made be- 
tween the two layers. We are particularly interested in 
the phase diagram in the cro-po plane, where po and (Tq 
denotes the total electron density and the imbalance pa- 
rameter (the normalized density difference between the 
two layers) , respectively. This is because phase transition 
points have experimentalljtl been observed most clearly 
in such a plane. We also present the phase diagram in 
the the po-AsAS plane, which will be useful to analyze ex- 
perimental data obtained in samples with different Asas 
but all other parameters unchanged. 

In Section II we review the Landau-site Hamiltonian 
governing bilayer QH systemsiSi. In Section HI we solve 
the eigenvalue equation for the ground state in the SU(4)- 
invariant limit. In so doing we identify the physical vari- 
ables within 15 generators of the group SU(4). We show 
that there are 9 independent physical variables, among 
which 8 variables describe real Goldstone modes on the 
SU(4)-invariant ground state. We then calculate the 
ground-state energy of the full system in the first order 
of perturbation. In Section IV we verify the existence 
of the spin phase, the canted phase and the ppin phase 
in the presence of the bias voltage. The ground state is 
obtained in an analytic form in each phase. In particular 
the spin component Sz and the imbalance parameter do 
are calculated as functions of the bias voltage for typical 
sample parameters. In Section V we present phase dia- 
grams in the Asas-Az plane, in the (jQ-pa plane and in 
the pq-Asas plane. In Section VI we show that the spin 
phase is not modified but the ppin phase is considerably 
modified by higher order perturbations. Thus, the spin- 
canted phase boundary is an exact result, as is consistent 
with the resulliS obtained from the exact diagonalization 
of a few-electron system. In Section VII we interpret the 
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experimental data^ based on the phase diagram in the 
CTQ-po plane. Section VIII is devoted to discussions. 



II. LANDAU-SITE HAMILTONIAN 

Electrons in a plane perform cyclotron motion under 
strong magnetic field B± and create Landau levels. In 
QH systems the electron position is specified solely by 
the guiding center X = {X, Y) subject to the noncom- 
mutative relation, [X, Y] = with the magnetic 

length £b = \/h/eBj_. It follows from this relation 
that each electron occupies an area 2tt£^ labelled by the 
Landau-site index. The total number of Landaus sites 
is = S/2'Ki% = N/v, where N and S are the total 
electron number and the area of the system, respectively. 
Electrons behave as if they were on lattice sites, among 
which Coulomb interactions operate. 

In the bilayer system an electron carries 4 different po- 
larizations associated with the ordinary spin (t,J,) and 
two layers (f,b). These four polarizations can be incor- 
porated in the isospin index /i =ft, fi, bt, b|. In this 
notation the second quantized electron field in the lowest 
Landau level (LLL) appear as 



m=0 



c^,{m){x\m), 



(2.1) 



where {x\m) is the one-body wave function for the Lan- 
dau site |m) in the LLL with m labeling the angular 
momentum. The operators Cp(m) and cl{n) satisfy the 
standard anticommutation relations. 

The electron field "0^ has four components, and the 
bilayer system possesses the underlying algebra SU(4). 
It has the subalgebra SUspin(2)®SUppin(2). We denote 
the 3 generators of the spin SU(2) algebra by t^p"^, and 
those of the pseudospin SU(2) generators by tPP'". There 
are remaining 9 generators in SU(4), which are given by 
r^P'^^T^'"". Their properties are summarized as 



^ppin^ppin ^^^^j ^ ^^^^^^ 



.ppin 



_spiii ppm _ ppin spin 



(2.2) 



where I stands for the 4x4 identity matrix. Their explicit 
form is given in Appendix IXI 

All the physical operators required for the description 
of the system are constructed as bilinear combinations of 
ip{x) and tp^x). They are 16 density operators 

p{x) ^'ip\x)^/j{x), 

Saix) =i^t(a,)^^pin^(a,)^ 
Paix) =^^j\x)tPP"^^j{x), 



(2.3) 



where Sa describes the total spin, 2Pz measures the 
electron-density difference between the two layers. The 
operator Rat transforms as a spin under SUspin(2) and 
as a pseudospin under SUppin(2). 

The total Hamiltonian consists of the Coulomb, Zee- 
man, tunneling and bias terms. The Coulomb interac- 
tion is decomposed into the SU(4)-invariant and SU(4)- 
noninvariant terms, 

H+ =i J d^xdhjp{x)V+{x - y)p(y). 

He =2 / d^xd?yP,{x)V-{x - y)P.{y), (2.4) 



where 



87re \ |a;| ^|a;|2 + ^2 



(2.5) 



with d being the layer separation. The tunneling and 
bias terms are summarized into the pseudo-Zeeman term. 
Combining the Zeeman and pseudo-Zeeman terms we 
have 

i?zpz = - jd^x (Az5^ + /^sAsPx + AbiasP^) , (2.6) 

with the Zeeman gap Az, the tunneling gap AgAS f^nd 

the bias term Abias = eVbias- 

Substituting the field expansion H2.1|l into we ob- 
tain Landau-site HamiltoniansiS., 

mnij 

He =4 E V^mjPAm,n)P,ii,j), (2.7) 

mnij 

where the Coulomb matrix element is 



y^mj - 2 ld'xd'y{m\x){x\n)V^ix-ymy){y\j), 



and 



(2., 



p(m,n) =E4("^)c<T(n), 

(T 

Sa{m,n) =iE4M(^r")..c.(n), 

(TT 

Pa{m,n) =iE4M(^aP^'")-c.(n). (2.9) 

(TT 

They satisfy the Woo (4) algebraili, which is the SU(4) 
extension of the Woo algebra. The Zeeman and pseudo- 
Zeeman terms read 

Hzpz = - '^[AzSz{n,n) + AsAsPx{n,n) + A'i.i^sPzin^n 

n 

(2.10) 

The total Hamiltonian is H ^ + iJ^ + Hzpz ■ 
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III. PHYSICAL DEGREES OF FREEDOM 

At v — 2 there are two electrons in one Landau site. 
We consider a creation operator of a pair of electrons at 
site n, 



(3.1) 



where is an antisymmetric complex matrix, g^i, 
— g^^. A homogeneous state is given by 



\g)=l[GHn)\0), 



(3.2) 



when g^^i, is independent of the site index n. The nor- 
malization of the state, (g|g) = 1, leads to 



Tr (ggt) ^ 2. 



(3.3) 



Since g^^, contains 6 independent complex parameters, 
there are 6 independent homogeneous states. They 
span the 6-dimensional irreducible representation of 
SU(4)(g)SU(4). 

It is hard to diagonalize the total Hamiltonian. Since 
we are interested in the regime where the SU(4)-invariant 
Coulomb term dominates all other interactions, we 
start with the ground state of the Hamiltonian , 

H+\g)^E+\g). (3.4) 

We include the SU(4)-noninvariant terms as small per- 
turbation. 
By requiring 



p{m,n)\g) = 1^5mn\g), 

the eigenvalue equation 13.4|l is satisfied withiS 



(3.5) 



and 



'X 



Here 



(3.7) 



(3.8) 



is the Coulomb energy unit. Note that the direct energy 
part, v'^ J2mi ^mmii^ Cancelled out by the background 
neutralizing chargei^. It is the ground state in the SU(4)- 
invariant limit of the system, which is the unperturbed 
system realized in the limits d ^ 0, Az — > 0, AgAs — > 
0. 



We introduce the expectation value of isospin opera- 
tors, 



§„=(g|^a(n,n)|g) = -Tr(r:P-5gt)^ 
5'„=(g|P„(n,n)|g) = lTr(rr'W), 



Jiab^{g\Rab{n,n)\g) ^^Tt{t'^ 



pin^ppin^^t 



(3.9) 



They are the total spin per one Landau site, and so on. 
We pay a special attention to the imbalance parameter 
(To, which is the normalized density deference between 
the two layers, |(To| < 1. At = 2 it is defined by 



_ Pf - Ph rp 

Pi + Ph 



(3.10) 



The relation (|3.9|) together with the normalization con- 
dition (|3.3|) yields 

99^ = 1^+1 + rPP- + T^rrr^'^ab) ■ 

(3.11) 

Here and hereafter the summation over repeated indices 
over the spin or the pseudospin is understood; for in- 
stance, §2 = 82 = J2a=xyz 

We study the condition (|3.5|) . It holds trivially for 
m = n. The condition p(m, n)|g) = for m ^ n yields 

<?a/35M-ci(m)4(m)ct (m)4 (n)|0) = 0, (3.12) 



which in turn leads to 



(3.13) 



where eai3tiv is the totally antisymmetric tensor. It is 
equivalent to 



1, 



(3.14) 



(3.6) in terms of physical variables, as we verify in Appendix 

o 

It is interesting to note that the magnitude of the 
SU(4) isospin, S"^ + 7^ is fixed by the eigenvalue 

equation (|3.4() . Namely it is a dynamical variable. In- 
deed, we have 



cos^ 261 



(3.15) 



for the state 

|g) = n [cos0Cf\(n)4(n) -Hsin0 4^(n)cj,^(n)] |0). 

n 

(3-16) 

In general, we can only derive a kinematical constraint 



+ y2 + Dl^ 1 



(3.17) 



from their definition H3.9|) : See ljB16|l in Appendix 
This is in sharp contrast to the v — 1 case, where the 
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magnitude of the SU(4) isospin, §^ + 3^^ + 31^ = 3/4, is 
fixed kinematically. 

Though there are 12 real parameters in the antisym- 
metric matrix g^^,, one of them is unphysical for repre- 
senting the overall phase, and another is fixed by the 
normalization condition H^.Hfl . Then one variable is fixed 
by the condition (|3.14|) on the ground state |g). As we 
shall see soon [see (|3.26f) ] and prove in Appendix lO there 
exists another unphysical variable which decouples on the 
ground state |g). Hence the parameter space character- 
izing the SU(4)-invariant ground state |g) contains 8 real 
independent variables. They are the 4 complex Gold- 
stone modes associated with a spontaneous breakdown 
of the SU(4) symmetry in the SU(4)-invariant limit. The 
number of the Goldstone modes agrees with our previous 
result^^'^^ obtained based on an effective theory with the 
use of composite bosons. 

We now include the SU(4)-noninvariant interactions 
as small perturbation. The first order perturbation is 
to diagonalize the full Hamiltonian H within this sub- 
space. Equivalently, the ground state is determined by 
minimizing the energy 

= £,N^ = (g|iJ|g) (3.18) 

within this parameter space. See Section IVIIII on this 
point. 

It is straightforward to show that 



(g|4(TO)c,.(n)|g) = J™„(5g^)^p 



(3.19) 



and 



14 (m)4(i)c^(j)c^(n)|g) 



for m 



i99^)Tcr{g9^)iyfi for m > 

{99'^)rii for TO > 

{99'^)y<y{99^)rii for < 

{99^)Ta{99^)yi^ for m<i,j>n 



,j <n 
J > n 
J <n 



(3.20) 



Thus 
1 



l-^T.^rnn.,{G\^lirnH{^)criJ)cAn)\G) 

= i99^)Ta{99^)u^e^ - (.9.9^)i'<T(gff^)rMex- (3-21) 

By using these formulas the ground-state energy per Lan- 
dau site is calculated as 

£g =e^ [Tr (rPP'"g4)]' - ^xTr (rr"55VPP'"55^) 
-e+Tr (.gsW) - ^AzTr (rr-g^) 



'2 AsAs 



Tr (rPP'"55^) 



-AbiasTr (rPP-5gt) , 



(3.22) 



AzSz- AsAsO'.- Abias?., (3.23) 



where is given by H3.7|l and 



(3.24) 



It is necessary to express 15 isospin components, §a, 
and 5iab in terms of independent variables. 

We may choose §a and Ta as 6 independent variables. 
We expect that Jiab are written in terms of these 6 vari- 
ables and 4 extra variables. Indeed, as demonstrated in 
the Appendix 1X1 we obtain 

■J<ab = ^ 7K7^ -X-PS 



§Q 



H — —■J<-PQ 



SQ 



Q 



Q 



7Q 



QS 



QaQb 

Q Q 



'-51 



QQ 



(3.25) 



with Q, = Sabchyc, § - |S|, 3^ = 1^1 and Q = |Q|. 
Here, 3ips, "Rpg, "^qs and ^qq are the 4 extra vari- 
ables. However, it can be proved that only 3 of them are 
independent, and they are parametrized as 



■31 



PS 



tJiQs -e'" ( -a + i§ ) , 



(3.26) 



Thus 3lab are expressed in terms of 9 variables; ^, A, 
Lu together with 8^, "Pa- One variable is unexpectedly 
unphysical, about which we explain in Appendix IBl 

We show that §x = §y ^ ^ij = for the ground 
state. First of all §a and 3iaz rotate as vectors under 
the SUspin(2) transformation. So, if we have any con- 
figuration with nonvanishing S^, and §y, we can perform 
an SUspin(2) rotation to increase S^, without affecting §^ 
and as far as possible until §x = §y = 0. This 

decreases the energy (|3.23|) . Similarly, performing an 
SUppin(2) rotation in the xy-pseudoplane, we can lower 
the energy via the tunneling term by increasing , with- 
out affecting and Dl^^, as far as possible until — 0. 

Substituting §x ~ §y = J'y = into Rab we come to 



q2 

^•^ab -72 + ^ 

ab ^ 



n + n: 



(3.27) 



Hence the energy H3.23|l yields 



(ex - 4) 



2£x§z 



ffcapJ'z 



AzSz 



2e 



^x ^x + 2ex rp2 ^ j)2 



J 



AsAsJ'a; " AbiasJ'z 



(3.28) 
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where we have set 



-cap 



4^0 -2£x- 



(3.29) 



Here ecapJ"^ is the capacitance energy per one Landau 
site. The capacitance parameter H3.29|l is different from 

v=l 
cap 



the ly — 1 QH systemifi, where ej^^^ = 4£^ — 4e^. Note 
that the energy formula (|3.23|) holds as it stands also at 
ly — 1. The difference arises because J2a-^az — -^z 
(pT^ at 1/ = 1 but it is given by (jT77| at 

If we minimize £g within the parameter space of the 
SU(4)-invariant ground state, we should impose the con- 
dition H3.14|l . However, it is instructive to minimize £g 
without requiring it. Namely, we only assume the kine- 
matical condition (|3.17|1 . which reads 



1 



^)§l + {J>l + J>l){l + e)+X^ <1- (3.30) 



Nevertheless, the equality is easily seen to hold on the 
minimum-energy state. Indeed, let us assume that it 
is realized with the inequality. However, this is self- 
contradictory since we can decrease the energy £g by in- 
in (|3.28|) as far as the inequality is obeyed; 



creasing 
note that ei 



The self-contradiction is resolved 



only if the equality holds in 13.30|l . 

(^1 + ^) §^ + in + n) (1 + e) + = i. (3.31) 

Namely, even if we minimize the energy £g without im- 
posing H3.14|l . we reproduce it. Consequently, the vari- 
ation approach reproduces the ground-state condition 
(|3.5fl within the class of test functions H3.2|l . 

We note that uj disappears both from H3.28|l and H3.30|l . 
and hence it is a zero-energy mode even in the SU(4)- 
noninvariant case. We shall see that it is related with 
the rotational invariance in the xy-plane: See (|4.3|l . 

We eliminate in by using (j3.31|) . 



?2 



1 



Az§. 



24 



2el 



(3.32) 



It is clear that we can decrease the energy by increas- 
ing §2 without affecting other terms in H3.32|l . This is 
achieved at by decreasing A^ until A^ = in H3.31fl , which 
yields 



We solve this as 



(3.33) 



1 
1 

7i+F 



a(3 



(3.34) 



in terms of two parameters \a\ ^ 1 and ^ 1. 
Substituting these into H3.32|l we obtain 



21 



Az^ r 5" AsAs 



^bias 

7i+l 



=a/3- 2e+ - 2e^. 



(3.35) 



To minimize this with respect to a, [3 and we write 
down the equations daE — dpE = d^E = 0. We rear- 
range them as 



ALs 4£x(A§-/3^Ai^s) 
1-/32 



Aov/r^ 



/3A, 



SAS 



Ao 



(3.36a) 



:, (3.36b) 



Az VI 
AsAS a 



(3.36c) 



where 



Ao 



^SAS 



Ai(l 



^) (1-/32). 



(3.37) 



The ground state is determined by these equations. The 
parameters a and (3 are solved out from H3.36a|l and 
(|3.36b|) in terms of the sample parameters [FIG^. Then, 
^ is given by H3.36c|) in terms of them. 



1 

0,8 
0,6 

0, ^ 
0,2 
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spin \ 


ppin 











1 

0,8 






' 
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\ ^/ / 


ppin 








0.2 










so 



100 



150 200 260 300 



FIG. 1: The spin component S2 and the imbalance parameter 
(70 = J'z are illustrated together with a as functions of Abias 
for typical sample parameters (d = 23nm, Asas = 6.7K and 
po = n X 10" /cm^). The spin phase (a — 0), the canted 
phase (0 < a < 1) and the ppin phase {a = 1) are realized 
for Abias < AS,,, Al\^ < Abia. < A^;;^, and A^^",, < Abias, 
respectively. There are 3 phases for n — lA but only 2 phases 
for n = 0.9. 
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IV. GROUND-STATE STRUCTURE 



and 



We discuss the ground-state structure as a function of 
parameters a and /3. Substituting (|3.36c|l into (|3.34|l we 
get 



Az 



S. (1 - a^) v/T^, 



^SAS 2 



^0 



(4.1) 



Recall that S^; = and Ty ^ 0. Using (|X26|) in 

(|3.25|) . we perform all the necessary substitutions and 
obtain 



A: 



SAS 



Ao 

A2 

Ao 

AsAS 



5iyy + iJi.,y = + Q,^ ^ ^ ^ \ - ^2 ^ 



VT^^x/T^^e*'", (4.2) 



while the rest components vanish, Jl^a — 0. We have 
expressed all the isospin components §a, 7 a and 3iab in 
terms of three variables a, (3 and uj. 

Using 2Sf = 8a + 3?a^ and 28^ = 8a - 3?az we find 



§1 = 




1 AsAS 






2 Ao 






1 AsAS 






2 Ao 



— cP-\J\ — cosw, 



smtj. 



(4.3) 



We see that u> describes the orientation of 8^ and Sj^ in 
the xy-plane. It is the zero-energy mode associated with 
the rotational invariance in the ccy-plane. 

Employing the formulae derived in Appendix El the 
matrix g can be reconstructed from 1)4.1(1 and 1(4. 2|l . For 
the sake of simplicity we take cj = ^ and come to 

9 = 7T5f + llal + 7!5l + + lldl + ihl. (4.4) 
where 



1 



_spin\ ppm 



9l =2^^-'""'^' 



gl ^^V2rf "rPP-, 



„s /'i _ spin\ ppm 



p _I sp 



(4.5) 



2^/2 

7o=0, 

2V2 



v/l-a2 1 



Ao 
Az 



7r 



2 0-^^ 

As 



7o^=-^^yW^, 



V2 A, 







ia ( AsASfl 



(4.6) 



Equivalently, the ground state is expressed as 



l{n)) =7fl§TH)+7ol§o(n))+7ll§i(n)) 

+ l^\Un))+ll\%{n))+ll\y^.{n)), (4.7) 



where 



I§t(" 

|8o(ri 
l§i(" 

|To(n 
|5'b(n 



==Cf\(n)c|,^(7i)|0), 
V2 

=Cf\(n)ct^(n)|0), 
=Cf\(n)Cf\(n)|0), 

=c|;^(n)c|;^(n)|0). 



|0), 



|0), 



(4.8) 



This reveals the microscopic structure of the ground 
state. 

All quantities are parametrized by a and (3. When we 
solve ((3.36all and ((3.36b(l for them, we find a < for 
Abias < A^'r^g, and a > 1 for Abias > A^Jg^^ with cer- 
tain values of Af^'j^^g and A{^|l^g: see (15.711 . Accordingly 
the ground-state energy ((3.35(1 is minimized by a = for 



Av 



< A^bias: and a 



1 for Abias > AP' 



bias ■ 



Then, 



the spin component 8z and the density imbalance 7z 
are calculated from 1(4.1(1 . In FIG^we illustrate and 
CTo = together with a as functions of Abias for typ- 
ical sample parameters; d = 23nm, Asas ~ 6.7K and 
Pq — n X lO^^/cm^ with n = 1.4 and = 0.9. 

First, when a = 0, it follows that 8^ = 1 and 7^ = 
since Ao = Az\/l — P^. Note that /? disappears from 
all formulas in ((4.6(1 . The spin phase is characterized 
by the fact that the isospin is fully polarized into the 
spin direction with 82 = 1 and all others being zero. The 
spins in both layers point to the positive z-axis due to the 
Zeeman effect. Substituting a — and Ao = Azy/l ~ fS"^ 
into ((4.6(1 we find the ground state to be 



igspi„>=n4(")4Twio)- 



(4.9) 
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It is interesting to notice that, even if the bias voltage is 
apphed, no charge transfer occurs between the two layers 
(cto = 0) as far as Abias < A^'^jj^g, where the system is in | 
the spin phase [FIGQ. 

Second, when a = 1, it follows that = and 7^ 
(actually CP^ + = 1). The ppin phase is characterized 
by the fact that the isospin is fully polarized into the 
pseudospin direction with 



(4.10) 



and all others being zero. Because Tz represents the den- 
sity difference between the two layers, /? is identified with 
the imbalance parameter co- Substituting a = 1 into 
(14. ()|) we obtain the ground state as 



gppin) = n I ^ (Vl + CToCf\(n) + VI - croci^in) 



X (VTT^Cf\(7i) + VT^c|,^(n)) } |0). 

(4.11) 

It is found that all electrons are in the front layer when 
ao — 1, and in the back layer when (Tq = —1. 

For intermediate values of a (0 < a < 1) none of the 
spin and pseudospin vanish, where we may control the 
density imbalance by applying a bias voltage as in the 
ppin phase. The ground state |gcant) is given by (|4.7|l 
with H4.6|) . All states except |§o('^)) contribute to form 
the ground state. This is so even in the balanced con- 
figuration (/3 = 0). It follows from H4.1(l and H4.8(l that, 
as the system goes away from the spin phase, the spin 
begin to cant and make antiferromagnetic correlations 
between the two layers. Hence, it is called the canted 
antiferromagnetic phasei^. 

We conclude that there are three phases in general; the 
spin phase, the canted phase and the ppin phase. They 
are as characterized by 



spin 


canted 


ppin 


= 1 




= 


J"^ = 




^2 = 1 



(4.12) 



The order parameters are and IP^ . 

Let us study the system with no bias voltage in more 
detail. The spin and pseudospin are 



A, 



SAS*- 



Ai(l 



(1 - a^) 



A: 



SAS 



1^2 + A| (1 



V^SAS" ^ 

where is easily obtained from (|3.36a(l with /? = 0, 



(4.13) 



^SAS 

4(2£x)2 



^Z 



A2 - A2 ■ 

^SAS ^Z 



(4.14) 
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canted 
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FIG. 2: The ground-state energies 



spin cppin 



and £i 



given in the spin, ppin and canted phases in the balanced 
configuration by taking d = 0.56£b. The horizontal axis is the 
tunneling gap Asas in the Coulomb unit Eq — /{ATvets)- 
A phase transition occurs along the heavy curve continuously 
from the spin phase to the canted phase, and then to the ppin 
phase, as Asas increases. The phase transition points Ag^s 
and Ag^g are given by 14.171 and 14.181 in the first order of 
perturbation. The energy is not modified but S^'^"^ is 

modified into fviC" by higher order quantum corrections. The 
point A stands for the would-be crossing of the spin and ppin 
energy levels in the first order of perturbation, with Aa given 
by 14.161 . The point B stands for the one when higher order 
quantum corrections are taken into account, with As given 
by 16.13L suggesting that the canted phase is considerably 
shrunk. 



The ground-state energy in each phase is given by 



^r" = -24-2ex-Az, 



24 



A2 A2 
^SAS I ^Z 



8e; 



8e; 



^^X^SAS 

A2 - A2 ' 
^SAS ^Z 



2e+ - Asas- 



(4.15) 



We have depicted them as a function of Asas in FIG El 
If the canted phase were ignored, the two levels £|p™ 
and fPP'" would cross at 



ASP 
^SAS 



2e 



X' 



(4.16) 



and a transition would occur suddenly from the spin 



phase (S 



1. 



0) to the ppin phase (§ 



0, 



T=l). Actually, a mixing of states occurs and lowers 
the ground state energy. As a result the level crossing 
turns into an level anticrossing [FIG [2], and the canted 
state emerges between the spin and ppin phases. As 
Asas increases, the phase transition occurs continuously 
from the spin phase to the canted phase, and then to the 
ppin phase. The phase transition points are 



from £|P" 



A SC _ 

^SAS — 



5^^"*, and 



A pc _ o - 
^SAS — ^^X 



4exAz 



A| + (2£x)^ 



(4.17) 



(4.18) 
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from £PP 



ccant 



These two formulas agree with the 
variational result due to MacDonald et alM; where the 
variational state has been chosen from a subset of the full 
set satisfying 



V. PHASE DIAGRAMS 

The phase diagram can be studied based on analytic 
formulas given in the previous section. We first present 
the phase diagram in the AgAS-^z plane for typical val- 
ues of Abias to compare our result with the standard 



ones 
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We search for the boundaries separating the 
canted phase from the spin and ppin phases in the sys- 
tem. These are extracted from H3.36a(l and Ij3.36b|l . The 
merit of our formalism is that analytic expressions are 
available to determine the phase boundaries even for im- 
balanced configurations. 

Along the spin-canted boundary we have a = 0. Sub- 
stituting this value into (|3.36bp we get 



/3 



A^ 



(5.1) 



v/l-/32 AsAS Az + 4ex ' 
We solve this for (3 and substitute it into (|3.36a|l to get 



A^AS = A| + 4exAz - 



Az + 4ex ■ 



(5.2) 



This determines the spin-canted boundary. 

Along the ppin-canted boundary we have a = 1. Sub- 
stituting this value into (|3.3t)a|l and Ij3.36bp we get 



iSAS 



2ec 



and 



A2 - 

^z — 



A 



bias 



2ec 



f3 



8s; 



(5.3) 



(5.4) 



after some manipulation. These two equations give a 
parametric representation of the ppin-canted boundary 
in terms of /3. 

In drawing the phase diagram in the Asas-Az plane 
we need to fix the layer separation d. Following the stan- 
dard literatur o^i^? , we have examined the case d = £b 
and presented the phase diagram in FIGOI Our results 
agree qualitatively with results obtained numerically by 
Brey et al.™ and MacDonald et alM- for imbalanced con- 
figurations. 

Our analytic formulas reveal some new features not 
reported in literature. We have found some peculiar be- 
haviors for the ppin-canted boundary as in FIG0I which 
occurs for 2e^ < 3e^, or d < 0.75£b- When the bias 
voltage is increased, a small ppin region appears in the 
vicinity of the origin at 



.(1) 



2ec 



(5.5) 




0.2 <!tA Ci.e 0.9 1 U 02 OA ^.a OA 1 
TLmrSlihg i^^^^ iLruirin^^ gup ■. 

FIG. 3: The phase diagram is given in the Asas-Az plane 
by changing the bias energy Abias- Here we have set d = Ib, 
and taken the Coulomb energy unit Eq for the tunneling gap 
(horizontal axis) and the Zeeman gap (vertical axis) . It is ob- 
served that both the canted and ppin phases are stabilarized 
even for Asas = by applying the bias voltage. The dotted 
curves in the top left panel represent the exact diagonalization 
result due to Schliemann et ali- for the 12 electron system. It 
is observed that the ppin-canted phase boundary is modified 
considerably by higher oder quantum corrections. 



This region is isolated from the basic ppin region by the 
canted phase. By increasing Abias, the isolated and basic 
ppin regions come closer, and they merge at 



.(2) 



3e^ 



(5.6) 



Then the isolated region disappears. 

These phase diagrams are, however, not so useful to 
analyze experimental data, since we need many samples 
with different d to realize, e.g., d = £b, but this is impos- 
sible. It is more interesting to see experimentally^ how 
the phase transition occurs by controlling the total den- 
sity po and also the imbalance parameter (Tq in a single 
sample with fixed values of d and Asas- Here we wish 
to describe new aspects of the phase diagram from this 
point of view. 

We start with the spin phase in the balanced configu- 
ration with Abias = 0. (The spin phase realizes provided 
the electron density is more than a certain critical value.) 
When the bias voltage is applied (Abias 7^ 0) the param- 
eter /? becomes nonzero according to (|5.1|l . However, 
no charge imbalance is induced between the two layers 
(ctq = 0) as far as a = 0. A charge imbalance occurs only 
above a certain critical value A^'r^^^ of the bias voltage as 
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FIG. 4: The phase diagram is given in the Asas-Az plane by 
changing the bias energy Abias- Here we have set d — OAEb, 
and taken the Coulomb energy unit Eq for the tunneling gap 
(horizontal axis) and the Zeeman gap (vertical axis). A small 
ppin region appears in the vicinity of the origin for d < 0.75£b- 
This region is isolated from the basic ppin region by the canted 
phase. This is shown in the top left panel. The upper part 
of each panel is the enlarged image of the vicinity of origin 
where the isolated ppin region appears. 



in FIG^ which is given by solving (|5.2II as 



\ bi; 



bias/ 



(Az+4ex)' 



^SAS- 



(5.7) 



This gives the spin-canted phase boundary in the Abias- 
po plane [FIGE]. For Abias > the system is driven 

into the canted phase with a ^ 0. As the bias voltage 
increases above a second critical value A^J^j^g, the system 
turns into the ppin phase with a = 1 as in FIG^ The 
critical value A^'^^^ is obtained by eliminating f3 in (|5.3(l 
and (|5.4|) . This gives the ppin-canted phase boundary in 
the Abias-po plane [FIG|S]. We remark that, as the elec- 
tron density decreases, the critical point A^'r^^^ decreases 
and eventually becomes zero so that the spin phase disap- 
pears at all: The critical density is po = 1.19 x 10^^/cm^^ 
in the case of FIG|5| 

To compare the experimental data^'^'^ it is more con- 
venient to present the phase diagram in the cro"Po plane 
[FIGini, which is constructed by using the relation be- 




BQ IW> 300 



FIG. 5: The phase diagram is given in the Abias-po plane for 
the sample with d = 23nm, where po — n x 10^^/cm~^. The 
two curves stand for the spin-canted phase boundary and the 
ppin-canted phase boundary in the first order of perturbation. 
The four large points are experimental data indicating phase 
transition points taken from Sawada et al?-. 




0.2 0.4 0.6 0.8 

Imbalance parameter oh 



FIG. 6: The phase diagram is given in the cro-po plane, where 
the sample parameters are the same as in FIG. |S] The spin 
phase is realized only in the balance configuration (ao = 0) 
along the heavy line. The solid curve represents the ppin- 
canted phase boundary in the first order of perturbation. The 
four large points are experimental data indicating phase tran- 
sition points taken from Sawada et ali. The dotted curve is 
a speculated ppin-canted phase boudary. The dotted hori- 
zontal line at n = 0.9 corresponds to the plateau-width curve 
indexed by n = 0.9 in FIGlHIb). 



tween the bias voltage and the imbalance parameter im- 
plied by and 

It is also useful to study the phase diagram in the po- 
AsAS plane, since it is not difficult to prepare samples 
with different Asas with all other parameters unchanged 
[FIGQ. This is constructed from iflTTjl and 
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FIG. 7: The phase diagram is given in the po-^SAS plane 
at the balance point, where the sample parameters are the 
same as in FIG. The horizontal axis denotes the total 
density po = n x 10^^ /cm^, while the vertical axis denotes 
the tunneling gap Asas in Kelvin. The dotted curve is a 
speculated ppin-canted phase boundary. The two large points 
are experimental data indicating phase transition points taken 
from Sawada et alii. 



VI. HIGHER ORDER CORRECTIONS 

We have explored the ground-state structure in the 
first order of perturbation theory. Our Hamiltonian con- 
sists of the SU(4)-invariant term and the SU(4)- 
noninvariant term i?non — + Hzpz- We are interested 
in the regime where :§> Hnon- Hence we have diag- 
onaUzed as in (|3.4(l and treat Hnon perturbatively. 
The eigenstates are degenerate with respect to , con- 
sisting of 6 states at each Landau site. The degeneracy 
is removed by the SU(4)-noninvariant terms i?non- We 
may construct 6 states \gi) satisfying 



(g,l^^l& 



(6.1) 



as follows. First we determine |go) by minimizing (g|_ff |g) 
with the use of a general state H4.7|l . We then construct 
a Fock space made of 5 states that are orthogonal to 
I go). Within this space we can minimize (g|iJ|g), which 
determines the lowest energy state |gi). In this way we 
may construct 6 states |go), ■ • ■ , Igs) with Eq < Ei < 
■ ■ ■ < E^, which satisfy Htj.lfl . The full Hamiltonian is 
diagonalized within the first order perturbation theory. 

We have constructed the phase boundaries between the 
spin phase, the canted phase and the ppin phase. Here we 
recall the exact-diagonalization result of a few-electron 
system showing that the boundary between the spin and 
canted phases is practically unmodified but the bound- 
ary between the canted and ppin phases is considerably 
modified from the mean- field result [FIG|3]. In order to 
explain this we discuss effects due to higher order per- 
turbations. 

We start with the spin phase. It is important that the 
fully spin polarized state (|4.9I) is an exact eigenstate of 
the total Hamiltonian H. This follows from 



PzihjMspin) = Px{hj)\gspin) = 



(6.2) 



and 

Sz{i, i)\gspin) = 1- 

We obtain F|gspi„) = ff'^gspin) with 



(6.3) 



^spin 



(24 



2e; 



A, 



The eigenvalue agrees with the mean-field value in H4.15|l . 
Hence the mean-field equations (I3.36a|l (|3.36c() are ex- 
act ones in the spin phase, implying that the spin-canted 
phase boundary is not affected by higher order perturba- 
tions. 

We proceed to discuss the ppin phase. For simplicity 
we consider the balanced point, where 

Szii,i)\gppin) = 1, ^';^(^*)|gppin> = 0- 

However, since 



(6.5) 



Pzihj)\ep 



7^ for i j, 



(6.6) 



though {gppin\Pz{i 7j)ISppin) — 0, the state |gppin) is not 
an eigenstate of the capacitance term H^. Hence the 
fully pseudospin polarized state H4.11|l is not an eigen- 
state of the total Hamiltonian H even in the absence of 
the tunneling interaction (Asas = 0). 

It is practically impossible to carry out the second or- 
der perturbation since almost no eigenstates of the un- 
perturbed Hamiltonian are known. We make a vari- 
ational analysis. We decompose the total Hamiltonian 
H into two pieces, Hq and Hi, where Hq is the maxi- 
mal part of H one of whose eigenstates is |g] 



ppin/ 



We 

do not give their explicit forms here since they are very 
complicated. We obtain 



^^olSppin) — '^^l^'^lffppii 



with 



gPPin 



A; 



SAS 



-capt-'o 



(6.7) 



(6.8) 



This agrees with the first-order perturbation result (|4.15|) 
at CTo — 0. Hence the variational analysis surely presents 
a higher order correction. According to the standard pro- 
cedure, we minimize the total energy with the variational 
state 



|gppin) = (l + Ai/l)|gppin>, 



(6.9) 



where A is the variational parameter. Then the energy 
reads 



cppiii _ op 



Kgppinl^^flgp 



'^F'"(gppinl^^l|gppi„> - (gppinl^^l-f^O-ffllgppin) 



(6.10) 

The calculation is straightforward though quite tedious. 
To perform various integrals explicitly we approximate 
as 



V+iq) = l.exp(^-l\q\^]El 



V-{q) = ^cxp(--|qhi?o 



(6.11) 
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with = eV(47re£B) 
result is given by 



which is vahd for d/£B ^ 1- The 



cppm 



gppii 



16i% 



1 

V2 



2 

^/3 



AsAs/-Eg 



(6.12) 

The correction is larger for larger layer separation d, and 
it is largest at the balanced point (ctq = 0). 

It is not easy to derive higher order corrections for the 
canted phase, since |gcant) is not an eigenstate of any 
simple Hamiltonian and furthermore it is a complicated 
state involving 5 states as in (|4.7() . Though we are unable 
to determine the improved ppin-canted phase boundary, 
we may make some arguments how it is modified. For 
simplicity we study the balanced point. We have shown 
that the would-be phase transition point is Ag'^g given by 
(|4.16l) in the ignorance of the canted phase. We examine 
how this point is modified by equating 1)6.4(1 and H6.12|l . 
Namely, from £^p"' = we find 



ASP 

^SAS 




(6.13) 

It is observed that this is much smaller than the first or- 
der perturbation result [FIGE]. Since the spin-canted 
boundary ((4.17(1 is not modified we expect that the 
canted phase is shrunk considerably. These features are 
precisely what are found in the exact diagonalization of 
a few electron system.^ [FIGE]. 



VII. EXPERIMENTAL STATUS 
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FIG. 8: The activation-energy data for po and the plateau- 
width data for ctq are taken from Sawada et al.^. In the upper 
figure (a), the point n — 0.9 is identified with the spin-canted 
phase transition point. In the lower figure (b), the minimum 
point in each curve indexed by n is identified with the ppin- 
canted phase transition point. The heavy curve traces these 
minimum points. 



Based on our theoretical results we wish to interpret 
the experimental data due to Sawada et el^^ yielding an 
unambiguous evidence for phase transitions. They have 
presented the activation-energy data for po [FIG|HIa)] 
and the plateau-width data for (Tq [FIG|HIb)]. In this 
section we set po = nx 10^^/cm^, and use n to represent 
the total density. 

We see a phase transition point (aQ,n) — (0,0.9) 
from the activation-energy data [FIGlH^a)], for which 
d/£B — 1-23 in the sample with d ~ 23nm. We give the 
phase diagram in AgAs-Az plane, upon which we plot 
the data point [FIG|H]. On one hand, the point is al- 
most on the spin-canted boundary. Furthermore, though 
the spin-canted boundary is exact, it is obtained in the 
ideal two-dimensional system and it will be modified in 
actual samples with finite quantum wells. Then it is rea- 
sonable that the point is slightly off the theoretical es- 
timation. On the other hand, as we have argued, the 
ppin-canted phase boundary is considerably modified so 
that the canted phase occupies a tiny domain in the phase 
diagram [FIGl21and FIG|3]. We cannot compare the ex- 
act diagonalization result^ with the data directly since it 
is available only at Vbias = and d/is = 1. Since the 



phase diagram for d/ls = 1-23 is similar to the one for 
d/is — 1, it would be allowed to extrapolate the exact 
diagonalization data from d/is = 1. We have plotted 
it on the same figure as indicated by the dotted curve 
[FIG|n]. Then the point is also near to the ppin-canted 
phase boundary. It is hard to decide on which phase 
boundary this point exists from this phase diagram. 

We answer this problem based on the phase diagram 
in the ao-po plane [FIGEJ. We focus on the horizontal 
line at n ~ 0.9. This corresponds to the curve indexed 
by n = 0.9 in the plateau- width data [FIGlUb)]. The 
curve takes a maximum at ctq = and a minimum at 
(To — 0.2 as the density imbalance is made. It is rea- 
sonable to identify the region ctq > 0.2 with the ppin 
phase. Then the region < ctq < 0.2 must be the canted 
phase. Note that the spin phase is realized only at o-q = 0. 
We conclude that the point (cto, n) — (0.2, 0.9) represents 
the ppin-canted phase transition point and that the point 
((To, n) ~ (0, 0.9) found in the activation-energy data rep- 
resents the spin-canted phase transition [FIG|S]. This in- 
terpretation is consistent with the fact that no indication 
of phase transion is found for 0.9 < n < 1.6 in the same 
data [FIGlHIa)]. 
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FIG. 9: The phase diagram is given in the Asas-Az plane 
for Abias = and d = 1.23£b- We have taken the Coulomb 
energy unit Eq for the tunneling gap (horizontal axis) and 
the Zeeman gap (vertical axis). The large point is the exper- 
imental data indicating a phase transition point taken from 
Sawada et ali. The two solid curves represent the phase 
boundaries in the first order of perturbation. The dotted 
curve stands for the ppin-canted phase boundary extrapolated 
from the exact diagonalization result obtained at d = Ib for 
the 12 electron system-. It is hard to decide on which phase 
boundary this point exists from this phase diagram. 

We continue to examine the plateau-width data for the 
sample with various density n [FIG|HIb)]. The data with 
n = 1.4 and n = 1.2 show similar behaviors as the one 
with n = 0.9. Namely, they possess the minimum points 
to be identified with the ppin-canted phase transitions 
around the imbalance parameter ctq ~ 0.3. It is observed 
that the ppin-canted transition point (Tq decreases as n 
decreases. It is not clear from the data when the ppin- 
canted phase transition disappears as n decreases. How- 
ever, it seems that the entire region belongs to the ppin 
phase a± n = 0.6. Let us tentatively regard the critical 
point exists near the point n = 0.6. In this way we have 
speculated the real ppin-canted boundary in the phase 
diagram in the (Jo-po plane [FIGE] and in the po-^SAS 
plane [FIG|7|. 



VIII. DISCUSSIONS 



and agree qualitatively with numerical results for nonzero 
bias voltages^ii^i. We have also argued how the phase di- 
agram is modified by higher order quantum corrections. 
We have shown that the spin-canted phased boundary 
is not modified but the ppin-canted phase boundary is 
considerably modified. It is necessary to develop a reli- 
able theory including higher order quantum corrections 
to determine the accurate ppin-canted boundary in vari- 
ous phase diagrams. 

Our new contribution is the phase diagrams in the uq- 
po plane as well as in the Abias-po plane. We have an- 
alyzed the relation between the imbalance parameter ao 
and the bias voltage Vbias- M v = \ the density imbal- 
ance occurs as soon as the bias voltage becomes nonzero. 
However, this is not the case at = 2 as in FIG^ where 
Abias = eT4ias- As far as 14ias < KTas density im- 
balance is not induced and the system is still in the spin 
phase. A charge imbalance occurs only above a certain 
critical value V^f^^ of the bias voltage, where the sys- 
tem is driven into the canted phase. As the bias voltage 
increases above a second critical value V^^^^s' system 
turns into the ppin phase. Taking these facts into account 
we have constructed the phase diagram in the tro-po plane 
as well as in the Abias-po plane. These phase diagrams 
will be useful to identify the three phases experimentally 
by using a single sample. We have also presented the 
phase diagram in the po-AsAS plane. It will be useful 
to identify the three phases experimentally by using sev- 
eral samples with different AgAS but all other parameters 
unchanged. 

We have interpreted experimental data due to Sawada 
et aliii^. We admit that our analysis is far from satis- 
factory. This is mainly because of the lack of qualified 
data. Recall that the primary concern of these experi- 
ments were to reveal the existence of the interlayer co- 
herent phase at v = 2^ which is identified with the ppin 
phase in the present terminology. In passing detailed 
experiments with several samples are urged to be per- 
formed in order to establish the ground-state structure 
in the v — 2 bilayer QH system. 



There are three phases, i.e., the spin phase, the canted 
phase and the ppin phase at u — 2. The ground state is 
the fully spin polarized state (|4.9|l in the spin phase, while 
it is the fully pseudospin polarized state (|4.11|) in the ppin 
phase. In the canted phase, however, it is not simply 
made of a certain combination of these two states, as 
postulated in a phenomenological bosonic spin model^^, 
but made of more states in a very complicated way. We 
have constructed the ground states explicitly in the first 
order of perturbation. 

The phase diagram has been studied numerically 
within the Hartree-Fock approximation in the standard 
literatur<4i^*2ii^. In this paper we have presented analytic 
formulas for it based on perturbation theory. Our results 
in the first order approximation reproduce precisely the 
Hartree-Fock variational resulti^ at the zero bias voltage. 
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APPENDIX A: PHYSICAL VARIABLES 

The basis of the algebra SU(4) is spanned by 15 Hermi- 
tian traceless matrices. For our purposes it is convenient 
to choose them as r^^'", tPp™ and t^p'"t^'"" given by 



spin 




ppm 




(Al) 



where I denotes the 2x2 identity matrix, and are the 
Pauh matrices. 

A basis in the space of 4 x 4 antisymmetric ma- 
trices, comprised of 6 matrices, can be chosen as 
T^P'"T^P'"rPP'" and tPp'^t^p'^tPp'". Hence, introducing 
four tree-dimensional vectors A, B, C and £), we may 
parameterize the matrix g as 



.spin ppin 



+ {Ay + zB,)(-zrPP'") + [Cy + iDy){-iTf^) 



= -i[{A + iBjT^'p''' + {C + iD)T 
We come to 



ppin spin ppin 



(A2) 



T-spin /I o 



(A3) 



(A4a) 
(A4b) 
(A4c) 

(A5a) 
(A5b) 
(A5c) 



on Sa, and 'Jiab- The constraint IjAScfl is verified based 
on the identity 

Sacd^bhe —Sab {^chSde — Sce^dh) — ^ah {Scb^de — SceSbd) 

+ Sae {5cb5dh - 5chhd) , (A6) 

while HA4a|) and IIA4b|) trivially follow from IjASa^ and 
(IX5bll . 

We construct Jiab satisfying ljA5ap - ljA5cp . To satisfy 
the constraint lASa^ we use the two normalized vectors 



2-a 'b 

Comparing this with (|3.11|) we get 

§a —£abcAl)Bc, 
•Pq ^abcCbDcj 
^ab =AaCb + BaDb, 

from which we derive the constraints 

§a^ab — 0, 
^abn = 0, 

^a-^b ^acd^bhe ^ch^de = 







§Q ' Q ' 



(A7) 



which are orthogonal to Sa, and orthogonal one to an- 
other. To satisfy the constraint (|A5b|l we use the two 
normalized vectors 



Q ' 



(A8) 



which are orthogonal to CPf,, and orthogonal one to an- 
other. Hence we are able to expand Hab in terms of 
4 tensors made out of these vectors with 4 coefficients 
!Kp5, 3ipQ, ^Qs and ^qq as in H3.25() in text. We now 
substitute (|3.25|l into the constraint (|A5c(l to find that 



PS-J^QQ 



(A9) 



Consequently three variables are independent among 
lKp5, 3?pQ, ^Qs and ^qq. They are parametrized as 
in ifX^ . 



APPENDIX B: UNPHYSICAL VARIABLES 

We have identified 9 independent physical variables in 
the 12 parameters of the antisymmetric matrix g. In 
this appendix we identify the 3 unphysical variables. In 
so doing we derive the kinematical condition (|3.17() . i.e., 

g2 _^y2 _^^2 ^ 1. 

The matrix g is expanded in terms of 4 three- 
dimensional vectors A, B, C and D as in ljA2p . We 
investigate how they are given in terms of the physical 
variables. For this purpose we reverse the constraints 
(|A4a|l and l|A4bp . The most general expressions read 



Aa 
Ba 
Ca 





- (S3')§a 




§Q 




- (S3')§a 




§Q 




- (S3')5'a 










7Q 



1 



9d + 



(Bl) 



where the parameters 9a, b and Ca,b are restricted by 



9aCb — Ca^b 
9dCc — CdOc 



(B2) 
(B3) 



We substitute IjBip into IIA4cl) , and compare the resulting 
equation with (|3.25|) . In this way we obtain new restric- 
tions, 





Ob' 




' Oc 




^PS 


_Ca 


<B . 




_0d _ 




■^QS _ 


' Oa 


Ob' 




' Cc' 




^PQ 


M 













(B4) 
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We solve these as 



We denote 



^ ^PsCb — ^Qs(^B a ^Qs(^A — ^PsCa 

8 ' § ' 

8 ' '^^^ 8 ■ 

(B5) 

It is trivial to check that ljB3|l is automatically satisfied 
due to l|A9|l and l|B2|l . We count the number of inde- 
pendent variables. There are 9 physical variables; §a, 
J'a, ^ps, ^PQ, 3?Qs and ^qq with one constraint IIA9|l . 
There are 3 extra variables; 9a, Ob, Ca and C.b with one 
constraint ljB2|) . We now show that they are the unphys- 
ical variables. 

A well-known unphysical variable is the overall phase 
of the ground state. First we identify it. We consider two 
vectors {9a, 9b) and (Ca, C-b)- We rotate the two vectors 
{9a, 9b) and {(a,Cb) by a single angle (j), which is the 
overall angle common to them. It turns out that the two 
vectors {9c, 9d) and {CcCd) rotate by the same angle 
due to IjBSjl . Then, {Aa,Ba) and {Ca,Da) rotate in the 
same way due to (jA5 b|l . Finally, ljA2|l implies that the 
matrix g acquires the angle 0, which is the overall phase 
of the ground state. 

Another well-known unphysical variable is associated 
with the normalization of the ground state. The normal- 
ization condition H3.3|l is equivalent to 



= 2. 



(B6) 



Substituting HB1() into this and using HB5() . we obtain 



- {9aCa + 9bCb) {^qs^ps + ^qq^pq)- (B7) 

We make a change of variables. Let 7 be the angle be- 
tween the vectors {9 a, 9b) and {C,a,C,b)- Omitting the 
overall angle we may write 



{9A,9B)={9cos'l,-9s\n^ 
(Ca,Cs) = (Ccos^,Csin^ 



Next, we introduce new variables x and y by 



(B8) 



Hps I ^2 
^QS + ^QQ 



QS 



QQ 



82 



^2 I fX)2 
Jipg JipQ 



'T)2 I ^2 



82 



(B9) 



which we use instead of 9"^ and C^. Now, there are three 
variables x, y and 7 with two constraints (jB2p and (|B7|I . 



82 (82 + ^2 + 3?2) 



or 



Here 



+ ^QS + ^Qq) + ^PS + ^PQ 



a 82 (82 + y2 + 3^2) 

l-a~ {^Rqs^ps + Xqq^pq) 



rp2 _ rp2 _ rp2 , rn2 , rn2 , rn2 

•'^ = •'^ab — •'^PS + -^PQ + -^QS + -'^QQ' 



(BIO) 



(Bll) 



(B12) 



where the last equality follows from H3.25|l . 

Using the above notations two constraints l|B2(l and 
(|B7|I are rearranged into 



COS7 



sin7 



V(l-a)(x2-y2)' 

2V82 
y^a{x^ - 2/2) 



(B13) 



These are well defined since > y^ and < a < 1, as 
follows from definitions (IB9II and (IBIOI). It follows that 



{x - 2)2 



4(§2+3)2_^3^2-) 



(1 — a)(a;2 — y"^) a{x'^ — y"^) 

from COS27 + sin27 = 1, or 



(B14) 



{ax - 2) V a (1 - a) t/2 ^ 4 (1 _ «) (1 _ §2 _ g)2 „ 3^2) _ 

(B15) 

We ask the question whether (|B15I) admits solutions for 
X and y which leads to positive values of 9^ and via 
(|B9|) . The solvability condition turns out to be 



(BI6) 



This is the condition H3.17|l for the magnitude of the 
isospin. 

If the condition is satisfied, we may solve x and y from 
((BT5|) . and obtain 9, C an d 7 f rom Ip9|) and ((BT3)| . Solu- 
tions are not unique. Ea. (|B15|) defines an ellipse. Moving 
along the ellipse the parameters 9, C, and 7 take different 
values, and therefore lead to different vectors A, B, C 
and D. However, all these vectors lead to one and the 
same set of physical fields 8a, and Jlab- This is the 
unphysical variable inherent to the v — 2 bilayer QH sys- 
tem, which is commented below (|3.26|) . This unphysical 
mode decouples on the ground state since the ellipse is 
shrunk to a point, as we see in the following appendix. 



APPENDIX C: GROUND-STATE CONDITION 

The eigenvalue equation 13.5|l in the SU(4)-invariant 
system leads to the condition (|3.13l) . or 



(CI) 
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On the other hand, the variational ground state in the 
full SU(4)-noninvariant system is found to satisfy H3.14|l . 
or 



31^ = 1. 



(C2) 



We verify the equivalence of these two conditions. 

For this purpose, we substitute l|A2|l into l|(yl|l . and use 
(IB6II. to obtain 



= + C = 1, AB = CD. 



(C3) 



On the other hand, examining this step of transformation 
we see that (l(]l|l follows from (|(]3ll . Namely, (|('l|l is 
equivalent to l|C3|l based on the formula (|A2(I . 

The derivation of (|C2p from HC3|) is easy. Using (|A4p 
we express + 3"^ + 01^ in terms oi A, B, C and D. 
We then use to derive 

The derivation of (|C3|) from l|C2p is more complicated. 
When the condition ljC2|l holds, the ellipse defined by 
(|B15p is shrunk to a point set given by ax = 2 and y = 0. 
In this case, 9 and C are determined by solving (jB9p . 



9 — ^/S^ + 3^p5+pQ, 



C=^§2 + D?|s+3?|q, (C4) 
and ljB13l) becomes 

cos 7 — \/l — a, sin 7 = ^/a, (C5) 
where a is given by IjBlOp . First we note that 



A^+D^^9\ + Q\ + 9l+Cl, (C6) 



to which we substitute (jB5p and HB8|I . We then use (|B7|I , 
(|C4|) and l|C5|l . In this way we prove + = 1 after 
some straightforward calculation. Similarly we can verify 
all equations in (|C3|I . 

As we have seen, the unphysical variables 9, C, and 7 
are fixed in terms of the physical variables on the ground 
state, where the condition (|C2|I is satisfied. Thus, every 
ground-state configuration has a unique matrix g^^. We 
substitute (jX^ . and into (EH) and fUsjl . 

Choosing w = 7r/4, for simplicity, we find 



(C7) 



and 





7 


cos 






2 




7 


sin 






2 



2 VA^ + A|(l-/32)' 



(C8) 



Using fCTjl and (|C8|l in ||B8)| we obtain 



9a = a = vT^, 

A, 



= -Cs = - -r^ yi^yT^. (C9) 
^0 

Now, it is straightforward to express A, B, C and D in 
terms of a and /3, and hence to reconstruct the matrix 
as in H4.4|l . 
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